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FIDELITY OF THE INTEGRATED FORCE METHOD SOLUTION 

Dale Hopkins, Gary Halford and Rula Coroneos 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Surya Patnaik 
Ohio Aerospace Institute 
Brook Park, Ohio 44142 

The theory of strain compatibility of the solid mechanics discipline was incomplete since St. 
Venant s strain formulation’ in 1876. We have addressed the compatibility conditions for both the 
continuum and the discrete system. 2 This has lead to the formulation of the Integrated Force Method for 
finite element analysis and the Completed Beltrami-MichelTs formulation in elasticity. 4 IFM is easily 
specialized to obtain the existing methods. 5 The reverse course cannot be followed. For example, the 
stiffness method cannot be specialized to obtain IFM because the stiffness method formulation is based 
only on the information contained in three quarters of the pie chart, see figure 1. The IFM variational 
functional yields the boundary compatibility conditions that were missed since 1876. 6 ’ 7 The compatibility 
limitation while blocking the growth of the primal method caused a bifurcation into the stiffness method 
and the redundant force method as sketched in figure 1. The compatibility compliance by the stiffness 
method— attempted through mesh refinement— may not be achieved for difficult problems. The classical 
force method with ‘cuts’ and redundancy is cumbersome. A Dual Integrated Force Method with 
displacements as the primal variable has also been formulated. IFM and IFMD yield identical solutions. A 
stiffness method code with a small programming effort can be improved to reap many IFM benefits when 
implemented with the IFMD elements. 

A modest finite element code (IFM/Analyzers) based on the IFM theory has been developed. This 
code uses numerical integration but does not use reduced integration or bubble function techniques. For a 
set of standard test problems the IFM results were compared with the stiffness method solutions and the 
MSC/Nastran code in a controlled environment. 7 ’ 8 For the problems IFM outperformed the existing 
methods, overshadowing the simplicity of the IFM elements. Superior IFM performance is attributed to 
simultaneous compliance of the equilibrium equation and compatibility condition. A method may be 
fallible but a utilization of the compatibility concept can never hurt. The additional information on 
compatibility condition has benefited analysis, 7 ’ 8 structural design^ and optimization. 111 

Mr. Bob Harder of MSC/Nastran organization wrote to us expressing reluctance to accept the high 
fidelity IFM solutions. 11 He wrote: 

“Authors ... make wrong statements about the accuracy of MSC/NASTRAN in the paper.” 

We re-solved all the examples again. The versions of MSC/Nastran software used were- 70 7 0 
Silicon graphics/IRIX64 6.5 (2001), 69.0 IBM/RISC System/6000 (1997), and MSC/Nastran 1982 
Section 2.2 (1982, Ref. 12). 


No numerical inaccuracy was detected in any of the seven tables in the paper. 5 The exercise 
reconfirmed IFM performance for static and frequency analysis. A report was sent to Mr. Harder. 13 We 
did not receive a reply. 
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Dr. Halford conducted a peer-review on IFM and the response received is given reference 14. 
Under his guidance a website (http://sdwww.grc.nasa.gov/patnaik) 15 is being developed. The site 
addresses solid mechanics problems of the strength of materials curriculum. A user can compare three 
methods (IFM, IFMD and the stiffness method) for internal force, reaction, displacement, and stress for 
modest problems that can be modeled by about 5000 equations. A numerical test-bed has been developed 
to validate the code. It contains the cantilever problem that was also solved earlier. 7,8 This example uses 
the RC0405 IFM element with five force unknowns that support linear variation of normal stress (c x = d 
+ c 2 y and a y = c 3 + c 4 x) and constant state of shear stress (x = c 5 ). The results obtained for this problem 
are given in the table and graphed in figure 2. The earlier version of MSC/Nastran exhibited a residual 
error of 8.6 percent in the tip displacement of the beam (bottom insert in figure 2). But the solution from 
the current version of the MSC/Nastran code is identical to the IFM results (top insert in figure 2) The tip 
displacement and the diagonal terms of the matrix for HMD, two versions of MSC/Nastran code and the 
standard stiffness method are given in the table for a two-element model. The current version of 
MSC/Nastran code has now produced numbers that are identical to the IFM/IFMD solution. 



Method 
of force 
(IFM) 


Indeterminate structures 




Current version of MSC/Nastran exactly matches 
IFM , • 

Q. ~ 
••= C 

? C 



s f 

® CD 

An earlier version exibited 8.6 % residual error 

E D. 
O .2 
Z -o 



Number of elements, n 


Figure 1.- Compatibility barrier prevented the development of IFM Figure 2,-Solutions for a cantilever beam 


Solution to a two-element model 



IFM/IFMD 

Stiffness Method 

I 

MSC/Nastran, 1990 

MSC/Nastran, 2001 
(Identical to IFM/IFMD) 

Tip displacement in 10 3 in. 

-1.1008 

(Timoshenko:- 1 . 152) 

0.2576 

1.0055 

-1.1008 


Normalized with 10 6 

' — J 

Normalized with respect to IFM/IFMD so] 

iution 

Diagonal coefficient of 
governing matrix 

1(6.1) (16.6) (3.0) (8.3) 
(3.0) (8.3) (6.1) (16.6) 

(1.25) (1.03) (1.27) (1.04) 
(1.27) (1.04) (1.25) (1.03) 

'(1.0) (1.02) (1.0) (1.01) 
(1.0) (1.01) (1.0) (1.02) 

fa-0) (1.0) (1.0) (1.0) 
(1.0) (1.0) (1.0) (1.0) 


Industry recognizes the importance of analysis but research on new methods is not central to its worry 
because it is profit driven. It indeed has built magnificent structures in the past from the Great Wall of 
China to the Leaning tower of Pisa without the benefit of analysis. Industry will earn profit in the future 
even without investing in analysis methods. Industrial prosperity is directly linked to the economy of a 
society but has a superficially loose connection to the academic knowledge. Academia and research 
institutions, which traditionally researched the analysis methods, are becoming eager users of codes 
developed by few software houses. Because profit dictates the survival of the software business, it is 
unable to indulge in basic research. This conflict-in-interest stalemate situation can create a vacuum that 
may be very expensive to reverse, especially with the retirement of the experienced researchers. It is 
important to recognize that a solution to a solid mechanics problem is incorrect when it is compatibility 
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non-compliant. The researchers in the field must adjust existing solutions and generate new ones that 
comply with the new conditions. The premier academic and research institutions should endeavor to 
increase the basic knowledge in the field to avoid saturation and arrest the declining solid mechanics 
discipline. The intent of this communication is to disseminate new knowledge on the compatibility 
concept to the structural mechanics community. It is not directed as a criticism to any individual or 
institution. J 
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APPENDIX A 


RESPONSE TO MSC/NASTRAN 
1.0 MSC/NASTRAN INQUIRY TO AUTHORS 

Date: Tue, 16 Sep 1997 17:07:20 -0700 
From: rlh@macsch.com (Bob Harder) 

To: dale.a.hopkins@lerc.nasa.gov 
Subject: Problems with Journal Article 

Cc : tom.cully@plano.macsch.com, ron.dyer@plano.macsch.com, Charley . wilson@plano . macsch . 
corn 

TO: Surya Patniak, Rula Coroneos and Dale Hopkins 
FROM: Bob Harder (MSC) 

DATE: August 28, 1997 

CC: Ron Dyer, Tom Cully, Charley Wilson 
SUBJ : Problems with Journal Article 

REF: DYNAMIC ANIMATION OF STRESS MODES VIA THE INTEGRATED FORCE METHOD OF STRUCTURAL 
ANALYSIS. By Surya N. Patniak, Rula M. Coroneos and Dale A. Hopkins. IJNME 
Vol. 40, pp. 2151-2168 (1997) 

The authors show how to animate stresses, but make wrong statements 
about the accuracy of MSC/NASTRAN in the paper. Such conclusions as 
"Only IFM/IFMD analysers produces accurate stresses. Stresses predicted 
by MSC/NASTRAN are too low" are not true. Here are some examples. 

1. The data in the tables is must be incorrect. How can the frequencies for their 
"stiffness" method agree precisely in every case in Tables III, IV, and V with 
MSC/NASTRAN? It looks like it may be a transcription error. Table III shows 
MSC/NASTRAN results for the vibration frequency of an isotropic cantilevered beam. 
Their models CB_Model_2 and CB_Model_3 should both become MSC/NASTRAN HEXA8 with 
three elements, so how is it possible that MSC/NASTRAN gives such different results 
for these models? The frequencies given for MSC/NASTRAN do not agree with those 
found by MSC. Table III (f requency/beam theory) should be corrected as follows: 


Elements 

Result 

in paper 

MSC/NASTRAN 



Model 

Mode 3 

Model 

Mode 3 

6 

HEXA8 

1.637 

1.640 

1.0065 

1.0419 

3 

HEXA8 

1.014 

0.995 

1.0221 

1.2724 

3 

HEXA8 

2 . 672 

1.862 



3 

HEXA20 

1.041 

1.091 

1.0322 

1.01 

6 

HEXA20 

not reported 

1.0137 

0 . 9860 


It is not clear that the frequency should converge to beam theory. There are some 
corrections to consider. Higher order beam theory has transverse shear flexibility 
and rotary inertia, both of which will predict lower frequencies. The fully 
constrained root (not beam theory compatible) gives an added constraint which will 
produce higher frequencies. 


NASA/TM— 2002-21 1286 


5 



t ‘ ^ Tai l le , VI1 shows maximum vonMises stresses for a cantilever beam with a tip shear 
load. While MSC can not reproduce their numbers, I agree that their model will not 
get good answers. The reason is that they are using poor boundary conditions. In 
order to agree with beam theory, you must let a solid element model breathe at the 
constrained end. Using their coordinate system, you can constrain the x- 
lsplacement, but the y- and z-displacements should be constrained only in an 
average sense. This was done at MSC using an RBE3 element. 


vm stress Beam Theory: 720.0 
vm stress MSC/NASTRAN reported: 355.9 
vm stress MSC/NASTRAN mod BC: 718.7 
vm stress IFM/IFMD analyser: 629.7 


I conclude that with proper boundary conditions the MSC/NASTRAN I-IEXA20 is bettor 
than the IFM/IFMD analyser. 


3. Their analysis of the static deflection results in Table VII is based upon 
simple (Euler) beam theory. A better target value is based upon beam theory 
including transverse shear flexibility (Timoshenko) 

deflection Euler Beam Theory: 
deflection IFM/IFMD analyser: 
deflection 6 MSC/N HEXA20, reported 
deflection Timoshenko Beam Theory: 
deflection 6 MSC/N HEXA20, mod BC: 

Using their boundary conditions, we agree that MSC/NASTRAN yields 2.25e-3, however 
with the correct boundary conditions MSC/NASTRAN is closer than the IFM/IFMD 
analyser . 


2 . 304e-3 
2 . 270e-3 
2 . 250e-3 
2 . 3165e-3 
2 . 315e-3 


Bob Harder 
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2.0 HOPKINS RESPONSE TO MSC/NASTRAN 


X-Info: IDE / NASA Lewis Research Center 
X-Sender: smdale@popserve.lerc.nasa.gov 
X-Mailer: Windows Eudora Pro Version 3.0.1 (32) 

Date: Thu, 20 Nov 1997 10:09:12 -0500 
To: rlh@macsch.com (Bob Harder) 

From: Dale A Hopkins <Dale. A. Hopkins @lerc.nasa. go v> 

Subject: Re: Problems with Journal Article 
Cc: tom.cully@plano.macsch.com, ron.dyer@plano.macsch.com, 
charley.wilson@plano.macsch.com 

At 05:07 PM 9/16/97 -0700, you (Bob Harder) wrote: 

TO: Surya Patniak, Rula Coroneos and Dale Hopkins 

FROM: Bob Harder (MSC) 

DATE: August 28, 1997 

CC: Ron Dyer, Tom Cully, Charley Wilson 

SUBJ: Problems with Journal Article 

REF: DYNAMIC ANIMATION OF STRESS MODES VIA THE INTEGRATED FORCE 
METHOD OF STRUCTURAL ANALYSIS. By Surya N. Patniak, Rula M. 

Coroneos and Dale A. Hopkins. IJNME Vol. 40, pp. 2151-2168 (1997) 

...[snipped] 

TO: Bob Harder, MacNeal-Schwendler Corp. 

FROM: Dale Hopkins, NASA/Lewis Research Center 

CC: Ron Dyer, Tom Cully, Charley Wilson 

RE: Addendum to Journal Article IJNME 40:2151-2169 (1997) 

Dear Mr. Harder: 

Thank you for your earlier e-mail correspondence regarding the above-referenced journal article. Attached to this e-mail is our 
rePv in the form of a Technical Addendum which we will also submit to the publisher. If you are unable to open the 
attachment, which is a Microsoft Word for Windows 95 (Version 7.0a) document, I would be happy to mail or FAX it to you. 

I hope you will give us the benefit of the doubt that it was not our intention to give a negative portrayal of the MSC/Nastran 
product. In retrospect, I wish we had included the information contained in the Technical Addendum because I do feel it 
provides a more complete picture of the state-of-the-art and the proper place for our modest effort from an historical 
perspective on the development and evolution of finite element methods. If you would like to discuss this further, please do 

not hesitate to contact me again by e-mail or by telephone at 216-433-3260. 

Sincerely, 


Dale Hopkins 


Dale A Hopkins Phone:216-433-3260 

NASA/Lewis Research Center FAX: 216-977-7051 

21 000 Brookpark Rd MS 49-8 

Cleveland OH 44135-3191 E-mail: Dale.A.Hopkins@lerc.nasa.gov 
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3.0 REPORT SUBMITTED TO MSC/NASTRAN 


“Dynamic Animation of Stress Modes via the Integrated Force Method of Structural Analysis 

IJNME 40:2151-2169 (1997) 


“ e lssue . s raiseb by Bob Harder o f MacNeal-Schwendler Corp. are specifically addressed 
Issue T ResuSn Table U 7 ' Typ0graphical error in Table 3; lssu ® 2 - Results in Table 3; 


Elements 


he elemente available in the IFM/Analyzers code are standard elements with full integration 
schemes, i.e. with neither of the two advanced features available in MSC/Nastran elements- 1) 
reduced shear integration; nor 2) reduced shear integration with bubble function. Therefore in 
the journal article the IFM results were compared with MSC/Nastran results obtained using ’ 
elements without the advanced features. 9 

We have now solved the problems with and without the advanced features of the MSC/Nastran 
elements. For the 8-node MSC/Nastran elements, the results are labeled as: 

8-node standard isoparametric (2x2x2 integration) 

8-node standard isoparametric (3x3x3 integration) 

8-node reduced shear (2x2x2 integration) 

8-node reduced shear (3x3x3 integration) 

8-node reduced shear with bubble function (2x2x2 integration) - MSC/Nastran default 
element 


1. 

2 . 

3. 

4. 

5. 


The 8-node elements available in IFM/Analyzers (for the IFM and stiffness 
standard isoparametric integration. These elements are comparable to the 
element labeled as 1 above. 


methods) use 2x2x2 
MSC/Nastran 


For the 20-node MSC/Nastran elements, the results are labeled as 


a. 20-node standard isoparametric (2x2x2 integration) 

b. 20-node standard isoparametric (3x3x3 integration) 

c. 20-node reduced shear (2x2x2 integration) 

d. 20-node reduced shear (3x3x3 integration) - MSC/Nastran default element 


IUd ?v3rt d 4 e nrrH ntS availabl f in IFM/Analyzers (for the IFM and stiffness methods) use 4x4x4 

th» mo™ d ^ pa .?T e "S lntegratlon ' respectively. These elements are comparable to 
the MSC/Nastran element labeled as b above. 


Boundary conditions 


eieSI ™ /Analyzers . code does not currently have provisions for interpolation constraint 
® ®" e " ts ' noy ,dr ™ I,| -P° ln ! constraints, we have introduced another boundary condition to 
a ow the model to breathe at the constrained end. This boundary condition, labeled as 
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clamped , is described as follows. For the model with 20-node elements, the displacement u 
(along the beam axis) is suppressed for all eight boundary nodes. The displacement w (along 
the load direction) is specified to ensure symmetry about the y-axis in the yz-plane, by 
suppressing this component at nodes 13 and 14 as shown in Figure la. This boundary 
condition allows the beam to breathe along the z-direction. Likewise, the displacement v (in the 
y-direction) is specified to ensure symmetry about the z-axis in the yz-plane, by suppressing this 
component at nodes 7 and 9 as shown in Figure la. This boundary condition allows the beam 
to breathe along the y-direction. 

For the model with 8-node elements, the clamped boundary has the displacement u suppressed 
for all four boundary nodes. The displacement w is suppressed at nodes 1 and 2 as shown in 
Figure Id. This boundary condition allows the beam to breathe along the z-direction. Likewise, 
the displacement v is suppressed at nodes 1 and 5 as shown in Figure Id. This boundary 
condition allows the beam to breathe along the y-direction. The boundary condition used in the 
journal article, labeled as “rigid”, suppresses all three displacement components for all nodes at 
the boundary (see Figures 1b and 1e). 

At the tip of the beam, the 10 lb. load is distributed among the nodes of the models with 20-node 
and 8-node elements, as given in Figures 1c and If, respectively. 

The beam theory results for the cantilever beam are provided, for reference, as follows: 

• Maximum normal stress in x-direction is 720 psi 

• Tip displacement from Euler beam theory is 2.304x1 0 -3 in. 

• Tip displacement from Timoshenko beam theory is 2.31 65x1 0 -3 in. 

• First frequency is 224.825 Hz 

• Third frequency is 1405.125 Hz 

Issue 1 - Typographical error in Table 3 

Although Table 3, as submitted to the journal and as published in a NASA Technical 
Memorandum (Ref. 38), was correct, the published journal article contained the following 
typographical error. The first column, fourth row (associated with CB_Model_2) should have 
stated “Six/HX20_90” instead of “Three/HX08_90”. We regret any confusion this may have 
caused. 

Issue 2 - Results in Table 3 

The frequency results for the problem with both three- and six-element models are given in 
expanded Table 3a and 3b. Both rigid and clamped boundary conditions are considered. The 
five 8-node elements of MSC/Nastran (1 . - 5.) and the four 20-node elements of MSC/Nastran 
(a. - d.), as described earlier, are used along with 8- and 20-node IFM/IFMD and stiffness 
elements of IFM/Analyzers. 
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• Eight-node elements! The MSC/Nastran results in Table 3 as published in the 
journal article are identical to the results obtained for MSC/Nastran 8-node element 1 
(standard isoparametric [2x2x2 integration]), which is comparable to the 
IFM/Analyzers 8-node elements. Please compare column 4, rows 1 , 2, 5, and 6 in the 
published article with the MSC/Nastran element 1 results, given in columns 12 5 
and 6 in the expanded Table 3b. 

• Twenty-node elements: The MSC/Nastran results in Table 3 as published in the 
journal article are identical to the results obtained for MSC/Nastran 20-node element b 
(standard isoparametric [3x3x3 integration]), which is comparable to the 
IFM/Analyzers 20-node elements. Please compare column 4, rows 3, 4, 7, and 8 in 
the published article with MSC/Nastran element b results, columns 1,2, 5 and 6 in 
the expanded Table 3a. 

• Frequencies: The IFM/IFMD frequency results are in agreement with the results 
obtained with MSC/Nastran advanced elements 5 (8-node reduced shear with bubble 
function [2x2x2 integration], and d (20-node reduced shear [3x3x3 integration]). 

Issue 3 - Results in Table 7 

The stress and displacement results for the problem with both three- and six-element models 
are given in an expanded Table 7. Both rigid and clamped boundary conditions are considered. 
The four 20-node elements of MSC/Nastran (a. - d.), as described earlier, are used alona with 
the 20-node elements of IFM/Analyzers. 

• Stress Results: The normal stress in the x-direction calculated using the 
MSC/Nastran element b (standard isoparametric [3x3x3 integration]), which is 
comparable to the IFM/Analyzers elements, converged to 396.01 psi. for the 6- 
element model with clamped boundary conditions. For the same conditions, the 
IFM/IFMD element produced a normal stress of 715.77 psi, while beam theory gives 
720 psi. The results with the default (improved) MSC/Nastran element (with reduced 
shear [3x3x3 integration]) and the IFM/IFMD element agree with beam theory results 
with minor deviations. 

• Displacement Results: For the 6-element model with clamped boundary conditions 
the displacement results agree with Timoshenko beam theory as follows: 

• 1 .0000 for Timoshenko beam theory 

• 0.9869 for IFM/IFMD 

• 0.9836 for MSC/Nastran element b 

• 1 .0006 for MSC/Nastran element d (default element) 

• 0.9994 for MSC/Nastran element d with the RE3 constraint (the element 

referenced in the correspondence from MSC) 

• Boundary Conditions: The improvement in displacement results between the rigid, 
clamped, and MSC’s RE3 constraint boundary conditions for the 6-element, 20-node 
per element, model are as follows: 
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2.31 5x1 O' 3 (RE3, MSC/Nastran, element d [?]) 
2.31 6x1 O' 3 (clamped, IFM/IFMD) 

2.270x1 O' 3 (rigid, IFM/IFMD) 
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elements of the force and displacement methods (for 8-node, element 1 ; for 20-node 
element b). 
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20-Node Element 


8-Node Element 



u = 0 at all 8 nodes 
v = 0 at nodes 9 & 17 
w = 0 at nodes 13 & 14 

(a) Clamped boundary condition 


z 



u = v = w = 0 at all 8 nodes 

(b) Rigid boundary condition 
z 

7 tl9 8 


3 LJ A 

Load of (-7/6) lb at nodes 3,4,7, and 8 
Load of 4 lb at nodes 15 & 16 
Load of 10/3 lb at nodes 11 & 19 

(c) Load distribution 



u = 0 at all 4 nodes 
v = 0 at nodes 1 & 5y 
w = 0 at nodes 1 & 2 

(d) Clamped boundary condition 


z 



u = v = w = 0 at all 4 nodes 

(e) Rigid boundary condition 


z 



Load of 2.5 lb at nodes 3, 4, 7 and 8 


(f) Load distribution 


Figure Al.— Eight-node and twenty-node CHEXA elements with boundary conditions and load distributions. 
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Expanded Table 3a (1st and 3rd mode frequency normalized with respect to beam theory solutions) 



Six twenty-node element model (6HX20) 

Three twenty-node element model (3HX20) 


Rigid b 

oundary 

Clamped boundary 

Rigid boundary 

Clamped boundary 

Methods 

1st mode 

3rd mode 

1st mode 

3rd mode 

1st mode 

3rd mode 

1st mode 

3rd mode 

IFM/IFMD 

1.007 

0.983 

0.993 

0.968 

1.025 

1.070 

1.004 

1.038 

Stiffness 

1.014 

0.995 

0.994 

0.974 

1.041 

1.091 

1.006 

1.063 

MSC/Nastran 







a 



* 

* 

* 

* 

* 

* 

b 

1.014 

0.995 

0.994 

0.974 

1.041 

1.091 

1.006 

1.063 

c 


* 

* 

* 

* 

* 

* 

* 

d 

1.013 

0.988 

0.993 

0.965 

1.031 

1.016 

0.993 

0.976 


Expanded Table 3b (1st and 3rd mode frequency normalized with respect to beam theory solutions) 



Six eight-node elen 

nent model (6HX08) 

Three eight-node element model (3HX08) 


Rigid bi 

oundary 

Clamped boundary 

Rigid boundary 

Clamped boundary 

Methods 

1st mode 

3rd mode 

1st mode 

3rd mode 

1st mode 

3rd mode 

1st mode 

3rd mode 

IFM/IFMD 

1.005 

1.044 

0.998 

1.036 

1.021 

1.274 

1.007 

1.264 

Stiffness 

1.637 

1.640 

1.627 

1.628 

2.672 

1.862 

2.661 

1.854 

MSC/Nastran 


— r 



1 

1.637 

1.640 

1.627 

1.628 

2.672 

1.862 

2.661 

1.854 

2 

1.637 

1.640 

1.627 

1.628 

2.672 

1.862 

2.661 

1.854 

3 

1.064 

1.103 

1.046 

1.088 

1.090 

1.358 

1.058 

1.337 

4 

1.637 

1.640 

1.627 

1.628 

2.672 

1.862 

2.661 

1.854 

5 

1.005 

1.044 

0.998 

1.036 

1.021 

1.274 

1.007 

1.264 
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Expanded Table 7 



Maximum VorMVIises stress in psi 

Stress ct x at root = 720 psi 

Tip Displacement 

Euler 2.304E-3 in.; Timoshenko 2.31 6E-3 in. 

Kigia d 

ounaary 

Clamped 

boundary 

Rigid boundary Clamped boundary 

Rigid boundary Clamped boundary 

Methods 

6HX20 

3HX20 

6HX20 

3HX20 

6HX20 

3HX20 

6HX20 

3HX20 

6HX20 

3HX20 

6HX20 

3HX20 

IFM/IFMD 

629.7 

573.1 

713.1 

677.5 

757.6 

721.4 

715.8 

679.4 

2.270E-03 

2.214E-03 

2.316E-03 

2.286E-03 

btiTtness 









2.246E-03 

2.1 61 E-03 

2.314E-03 

2.279E-03 

MSC/Nastran 










a 

L 

1 .0E+06t 

1 .9E+05f 

1.1E+06t 

2.4E+05t 

799.0 

820.8 

719.7 

721.0 

ri.405E-02 

4.583E-03 

1.45 IE-02 

5.289E-03 

b 

355.2 

289.2 

399.4 

366.3 

404.3 

381.3 

396.0 

362.3 

2.246E-03 

2.161 E-03 

2.314E-03 

2.279E-03 

c 


5.7E+17t 

7.5E+17t 

2.8E+17f 

* 

3.4E+16t 

2.5E+16t 

-5.0E+16t 

* 

^350E+09 

3.749E+09 

2.892E+09 

d 

642.7 

520.9 

725.2 

722.8 

774.9 

778.6 

719.1 

719.6 

2.251 E-03 

2.192E-03 

2.318E-03 

2.318E-03 

Failed 






tHigh values 




















APPENDIX- 2 

External Reviews on Finite Element Method for the 21 st Century 

A specific question was raised: 

• ‘Should we invest in an IFM based Finite Element Method Project? or 

• You believe the stiffness method to be adequate even for the 21 st century. 

Excerpts from Reviews 


John Meek 
Australia 

Lucien Schmit 
USA 

John 

Przemieniecki 

USA 

George Thierauf 
Germany 

Hiro Miura 
USA 

“Stiffness 

“.. makes a 

“ IFM should be 

“Fully agree with 

“ ..intrigued 

method has 

persuasive case 

developed .. just 

.. advantages of 

by .. new 

reached the 

for further 

as NASTRAN 

IFM” 

compatibility 

plateau of 

development of 

was developed 

“Support revival 

equations..” 

diminishing 

IFM .." 

for 

of ‘natural’ way of 

return”. 


Displacement 

Computing 

“Develop 

“Rough 

“..Displacement 
methods to 

method..” 

stress..” 

IFM/IFMD.. 
Important to 

(stiffness 

have difficulty in 

“ accuracy 

“Plasticity: IFM 

solve 

method) .. 

predicting 

arguments in 

possesses 

1 million 

results taken as 

stresses 

favor of IFM are 

immense 

degrees of 

gospel truth - 
with religious 

accurately..” 

well known..” 

I advantages..” 
“Optimization 

freedom .. 
problem..” 

belief, which is 

“..invest in 

“Bottom line .. 

and 

dangerous”. 

IFM..” 

further research 
.. should be 

decomposition 
offers same 

“ ..New 
frontier in 

“ .. full and 

“..deserves high 

initiated by 

advantage as 

solid 

enthusiastic 
support., for 
project..” 

praise and 
respect..” 

NASA.” 

plasticity.” 

mechanics..” 

“..NASA to 
support 
prototype 
code 

development. 

» 


Seven reviews are enclosed. 
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National Aeronautics and 
Space Administration 

John H. Glenn Research Center 
Lewis Field 

Cleveland, OH 44135-3191 



Reply to Attn of: 


14* June 1999 


Re: Finite element method for the 21 st century. 


Dear—, 


As we know, only approximate stress solution can be obtained for the vast majority of structural 
mechanics problems. Thus the founding the fathers of this discipline preferred calculating stress 
by two competing methods; then ascertaining accuracy by comparing the two solutions. This 
simple reasoning is not currently applicable because only the displacement (or stiffness) method 
solution is available. Traditionally, accurate resolution of a stress state was not intrinsic to the 
stiffness method. Despite progress stress accuracy by the stiffness method cannot always be 
guaranteed. 

In this context we wanted to examine the potential of the integrated Force Method to 
become a possible alternate formulation to the stiffness method. Because of your comprehension 
of this field, we would like to receive your response to the question: Should we invest in an IFM 
based finite element method project, or do you believe the stiffness method to be adequate even 
for the 21 st century? 

We request your kind opinion on this matter. We look forward to hear from you at the 
earliest possible time. Let me thank you in advance for your time and consideration. 

Sincerely, 


Gary R. Halford, Ph.D. 

Senior Technologist 

M.S. 49-7, NASA Glenn Research Center 
21000 Brookpark Road, Cleveland, Ohio 44135 
Phone: (216) 433-3265 
E-mail: Garv.R.Halford@qrc.nasa.aov 


Enclosure: IFM report. 
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UNIVERSITY OF QUEENSLAND 

Civil Engineering Depar tment 

Professor J.L Meek B.E. B.S. M.S. PHD. 

President: Numerical Analysis Research Coorporation Pty. Ltd. 
Tel: (617) 3365 1569, or 3378 8423 
Fax: (617) 3365 4161, or 3378 8568 
Email: e2meek@brolga.cc.uq.oz.au 


St. Lucia 4072 

Queensland 

Australia 


May 18, 1999 


ATTENTION: 

Dr. Gary R. Halford 

Senior Technologist Your reference: Letter from Dr. Surya Patnaik 

M.S. 49-8 NASA Glenn Research Center 

21000 Brookpark Road 

Cleveland, Ohio 44135 

U.S.A 


Dear Dr. Halford 


Reintegrated Force Method 

Today I received a very interesting letter and document from Dr. Paraiak on the 
Integrated Force Method Analyser. The observation that the stiffness method 
has reached a plateau of diminishing returns is self evident to any who can look at 
the journal articles in an objective manner. The danger of the present situation 
is that results from rather rough finite element analyses are being taken to be 
gospel truth with a sense of almost religious belief. Thus if the integrated force 
method offers the possibility of a better description of the stress fields within a 
continuum then the method is worthy of development. Thus I give my full and 
enthusiastic support to the proposed research project and would like to be kept 
informed of its progress. Thank you for the opportunity for reviewing the project. 
With my very best wishes for its success, and with my very best regards, 


Yours sincerely 



John Meek 
Professor Emeritus Civil Engineering 
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L.A. Schmit 
1/4 


545 Third Ave., S. 

Edmonds, WA 98020 

May 25, 1999 

Dr. Gary R. Halford 
Senior Technologist 

NASA Glenn Research Center (M.S. 49-7) 

21000 Brookpark Road 
Cleveland, Ohio 44135 

Dear Dr. Halford, 

I am writing in response to your letter of May 5, 1999 regarding development of 
an industrial strength IFM analyzer”. While I have not done any significant finite 
element research for about 20 years, I find that the Presentation you sent me (entitled 
Finite Element Method for the 21 s1 Century”) makes a persuasive case for further 
development of the IFM in a finite element context. 

Historically force methods were favored in Civil and Mechanical engineering 
applications involving very stiff structures. Accurate stress information was needed and 
displacements were often so small that they were often not even calculated. Also the 
classical redundant force method was an easy to visualize natural extension of the force 
method as applied to statically determinate structures. 

In an aerospace context the stiffness method has long been favored because the 
structures to be analyzed are usually very flexible and they often involve finite elastic 
deflections, dynamic response, and elastic instability considerations. The displacement 
method (stiffness) of structural analysis has achieved a dominant position in aerospace 
structural analysis because is it well matched to the tasks that aerospace structural 
engineers generally face. 

However, it is important to recognize that the commonly used displacement 
methods dO-haye difficulty in predic ting stresses accurately . In my opinion, aerospace 
structural analysis problems where stress accuracy is of paramount importance, can 
benefit from application of the IFM. It is hard to get around the fact that in the commonly 
used displacement methods, stresses must be determined by taking first or second 
derivatives of approximate displacements. These differentiations amplify the errors that 
exist in the approximate displacement solution. 

In looking at various alternative methods in structural mechanics ( see copy of 
pl4 from the Presentation you sent me ) it is useful to ask which equations of elasticity 
end up being satisfied approximately when a discretized solution is sought. For example, 
in the displacement method the equilibrium equations of elasticity are satisfied 
approximately. This follows from the fact that the equilibrium equations are satisfied 
exactly” at only a finite number of mesh points (for purposes of this discussion I assume 
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L.A. Schmit 
2/4 


lhal ihe system displacements are truly single valued at all points). In the classical force 
method using stress functions (method 2 on p. 14, copy attached) it is the equations of 
compatibility that are satisfied approximately, because in a discrete formulation they are 
satisfied exactly at only a finite number of mesh points. 

Let us now ask which equations of elasticity end up being satisfied approximately 
when a discrete solution is sought based on the CBMF (IFM). The answer appears to be 
that in a discrete formulation the equilibrium equations and the compatibility conditions 
(the EE and the CC ) are satisfied “exactly” at only a finite number of mesh points, 
therefore both the EE and the CC are satisfied approximately. Since the complete 
Beltrami-Micheil Formulation is now available ( see attached copy of p. 13 from the 
Presentation ) as a set of partial differential equations expressed exclusively in terms of 
stresses, it would be interesting to seek a finite difference solution of the stress values at 
each mesh point. When this is done it would be best to use a test problem for which a 
closed form solution is available. 


In any event the bottom line is that a discrete CBMF solution generates stress 
valuesf at a set of mesh points) directly, therefore the errors associated with 
differentiation of approximate displacement distributions are avoided Furthermore, if 
strains or displacements are needed they can be obtained by back calculation (stresses -> 
strains -> displacements). The resulting strains and displacements will correspond exactly 
to the approximate stresses obtained via the discrete CBMF solution. 

The IFM appears to be particularly well suited to stress analysis problems 
involving the introduction of concentrated loads, composite material internal and edge 
stress evaluation, and fatigue evaluations where stress accuracy is important In the 
overall scheme of things I do not believe that the IFM will replace the current widely 
used displacement methods, but the IFM does seem to offer an essential n*w supplement 
to the currently available collection of tools for aerospace structural analytic 

In summary then, my answer to your basic question is YES, you should invest in 
an IFM finite element code project. However, it would be prudent to view the IFM finite 
element code as a high stress accuracy complement to existing displacement method 
programs, not a replacement for them. In closing I want to say that Dr. Patnaik deserves 
high praise and respect for successfully completing the Beltrami-Micheil formulation in 
elasticity. This accomplishment has academic significance independent of any 
subsequently developed finite element analysis capability. 

I hope that you will find my observations and comments useful in your 
deliberations. 


Sincerely, 

-ucien A. Schmit, Jr. 
Professor Emeritus, UCLA 
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Completed Beltrami-Michell Formulation in Elasticity 



CD-99-78688 


L.A. Schmit 
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Methods of Structural Mechanics and 
Associated Variational Functionals 


Method 


Primary Variables 


Method Elasticity 
Number 


Completed Beltrami-Michell 
Formulation (CBMF) 

Airy Formulation 
Navier Formulation (NF) 
Hybrid Method (HF) 

Total Formulation (TF) 


Structures 

Elasticity 

Structures 

Variational 

Functional 

Integrated Force 
Method (IFM) 

Stresses 

Forces 

IFM Variational 
Functional 

Redundant Force 
Method 

Stress Function 

Redundants 

Complementary 

Energy 

Stiffness Method 
(DM) 

Displacements 

Deflections 

Potential 

Energy 

Reissner Method 
(RM) 

Stress & 
Displacements 

Forces & 
Deflections 

Reissner 

Washizu Method 
(WM) 

Stress, Strains & 
Displacements 

Forces, Deformations & 
Deflections 

Washizu 

Functional 


• Methods (2 to 5) have been derived from IFM equations 



Dr. Gary R. Halford 
Senior Technologist 

MS 49-7, NASA Glenn Research Center 
21000 Broadpark Road 
Cleveland, OH 44135 

Dear Dr. Halford: 


22 June 1999 


I apologize for not responding earlier to the May 5, 1999 letter. I was away on a 
six week trip in Europe (England, Switzerland, and Poland where I received an honorary 
doctorate from the Warsaw University of Technology). I believe the proposal submitted 
by Dr. Patnaik merits support from NASA. Although the finite element methods 
originated from the Force Method, the competing Displacement Method quickly became 
the method of choice because of its simplicity. Today, because of the enormous 
computing capabilities available to us, there is no reason why the Force Method, or the 
generalized approach of the Integrated Force Method (EFM) should not be developed 
further and made available for structural analysis, just as NASTRAN was developed for 
the Displacement Method. The accuracy arguments in favor of IFM, presented by Dr. 
Patnaik are well known. Thus my bottom line is that further research into the Force 
Method should be initiated by NASA. 


You may find it interesting that in 1997, when I was a Visiting Professor with the 
University of Southern Florida (USF), I submitted a proposal dated 2 June 1997 to NASA 
Langley Research Center for the development of a general theory of the Force Method of 
analysis for aeronautical and space structures based on a new integrated formulation of 
equations of equilibrium and compatibility. The formulation of the governing equations 
was different from that of Dr. Patnaik, but the final results were the same. The proposal 
included a complete theoretical analysis and a pilot computer program, as a demonstrator 
The total cost of under S100K was to be shared between NASA (70%) and USF (30%). 
The proposal was turned down. It was argued that the state-of-the art of the commerciai 

fimte element analysis is such that the development of an alternative method could not be 
justified. 

Sincerely yours 



J. S. Przemieniecki. PhD, DSc 
President 

Astra Technologies, Inc. 


510 Pennyroyal Place, Venice, FL 34293-7233, Tei(94i)-496-9662, FAX <941)-496-9662 
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Prof, Dr.-lng. G. Thierauf 
Baumechanik - Baustatik 


feQ.f, Dr.-lng. G, Thie rauf -Universttat GH Essen - Postfach 1 0 37 fia . 4B rr? £< «»<■> 


Dr. Gary R. Halford 

M.S. 49-7, NASA Glenn Research Center 
21000 Brookpark Road 

USA-Cieveiand, OH 44135-3191 



Universitat 

Gesamthochschule 

Essen 


Fachbereich 10 Bauwesen 

UniversitatsstraBe 1 5 - D 45 1 1 7 Essen 


Vermittlung: 

Direktwahl; 

Sekretarist: 

Telefax: 

(intern,: 

E-mail: 

Raum: 


( 0201 ) 183 -0 
(0201)183 - 26 70 
(0201) 183 - 26 72 
(0201) 183 - 26 75 
0049/201 183-2675) 
th<8>bauwesen,un»~essen,de 
gThierauf@uni-essen.de 
VI 5 S03 D23 


Essen, 19.09.01 TH/pa 


Dear Colleague, 

some reason, which I cannot trace exactly (it might well be my fault} I received your letter 
dated June 3rd only today together with a second letter from Gary R. Halford. 

The question which you arise has directed much of my research for almost 10 years and I came 
up with the conclusion that we missed the train perhaps more than 30 years ago when the direct 
stiffness method started to dominate the force method. The very few researchers who went on 
working in the field of the force method and who contributed to the integrated force method 
(IFM) often had problems to justify their research. I fully agree with you with respect to the 
advantages of IFM and I would certainly support any attempt to find a revival of the .natural* 
way of computing stresses in structures. 

From my experience I would like to mention two other most favourable aspects: 

- Plasticity 

The IFM possesses an immense advantage when dealing with all variants of 
plasticity. In this field, where general non-unique functional relations between stress 
and strain exist, the primal formulation does not only offer a one order increase in 
accuracy but also a tremendous saving in computational work, 

- Optimization and decomposition 

Early attempts to exploid the advantages of the IFM date back to the early sixties. As 
most of our constraints in structural optimization are governed by stresses, a direct 
stress-based formulation as in IFM offers the same advantages as in plasticity. 

In addition, the IFM offers the possibility to generate an optimal vector basis for any iterative 
design. We have experimented with the basis generated by the compact eigenstresses and found 
out that we could overcome almost all the problems for decomposition in structural optimization. 

Unfortunately our proposals for funded research found little resonance in the German Research 
Association - the experts were convinced that the diplacement formulation was superior. So I do 
hope that you will be more successful with your attempt and I can assure you that I would be 
prepared to support it as much as possible. 

Sincerely, 


Georg Thierauf 

H:\Workarea\IFM_REVSEWS\Thierauf-Germany-1 1 .doc 
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National Aeronautics and 
Space Administration 

Ames Research Center 
Moffett Field, California 94035 


H. Miura 
1/3 



Reply to Attn, of Mail Stop 237-11 


Dr. Gary R. Halford 
Senior Technologist 
MS 49-7 

NASA Glenn Research Center 
21000 Brookpark, Cleveland, Ohio 44135 


Hirokazu (Hiro) Miura 
Aerospace Engineer 
Mail Stop 237-1 1 Systems Analysis 
NASA Ames Research Center 
Moffett Field, CA 94-35-1000 
Tel: (650)604-5888 
hmiura@mail.arc.nasa.gov 


June 11, 1999 


Dear Dr. Halford: 

In his letter dated May 5th, 1999, Dr. Surya Patnaik requested me to convey my 
opinion on the merits of the Integrated Force Method as the subject of the future 
NASA research. I have been aware of IFM for some time and had a nice 
conversation with Dr. Patnaik on this subject when he visited Ames about a year ago 
for the CAS meeting. However, I have never studied technical details of IFM or 
IFMD in the past. This time, I would like to study this subject more meticulously, 
but my out of country trip during May 15 through 23 and an intensive proposal 
preparation activity from May 28 through June 4th did not leave much time forme to 
do anything else. Since the deadline requested by Dr. Patnaik has passed already, I 
would summarize my impression on this subject in this letter. 

I was very much intrigued by the derivation of a new compatibility equations to 
complete the Beltrami-Mitchell formulation based on the new variational ftmction. I 
think it important to study the validity and implications of the content in the table 
summarizing five basic methods of structural mechanics from all the possible angles. 
Based on this table and the descriptions in the following articles, I interpreted that the 
contributions of IMF/IMFD in more detailed level may be identified as introduction 
of new compatibility equations in elasticity and the concept of working with 
"deformation variables" in structures. 
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I believe that Dr. Patnaik’s findings warrant systematic studies and developments to 
clarify a number of issues associated with IMF. Especially, basic theoretical aspects 
of IMF/IMFD in solid mechanics should be reviewed thoroughly, extended as 
required and documented in a self contained literature. On application side, we all 
know that commercial finite element structural analysis codes based on displacement 
methods have contributed a great deal to expand the horizon of structural analyses of 
structures of practical scale and complexity. Since these codes are now used as 
indispensable tools in the routine production environment, it will be difficult to 
change the way of doing business, unless very distinct advantages are found in the 
new approaches. However, there still exist such structures that cannot be analyzed 
satisfactorily or that are very difficult to create adequate analysis models within 
reasonable amount of effort. 

It will be very appropriate for NASA to provide support for the basic theoretical 
study and comprehensive review of IMF/IMFD. Also, it will be appropriate for 
NASA to support prototype code development to the extent that performance of 
analysis program based on IMF/IMFD can be measured quantitatively on the modem 
computational infrastructure, to make fair assessments compared to the capabilities 
and performance of the analysis codes based on displacement methods. 

In the following, I will try to list the subjects I thought of while reading the materials 
provided by Dr. Patnaik. 

1 . It will be helpful to formulate the weak form for finite element analysis based on 
the variational formulation of n t and to work out equivalence of strong and weak 
forms. 

2. Explore the possibility of expanding IMF/IMFD to boundary element analysis of 
solids. Granted that I do not have any background in BE, but I am aware that the 
surface stress of solids obtained by straight-forward BE formulation is far from 
accurate, i.e. unusable in practice. On the other hand, some of the complex solid 
structure can only be modeled for analysis by BE codes based on the surface 
description available from CAD data. Can IMF/IMFD extend any help? 

3. It is necessary to study basic characteristics and performances of the finite 
elements developed by IMF/IMFD. Especially, it will be very useful if accurate 
stresses can be obtained at the locations where stress gradients are large without 
using extremely fine mesh. Also, it is equally important to study the 
performance/robustness against element shape distortion from the ideal geometry. 

4. In practice, a method that cannot be extended to deal with dynamics and problems 
in nonlinear regime may not be able to enjoy widespread acceptance. The results 
on natural frequencies of the turboprop blade reported in NASA TM 4729 are 
rather bizarre and such items must be studied very much in detail. 

5. Design optimization capability is an increasingly important feature required for 
any structural analysis programs. In this context, it is important to perform 
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research on feasibility of efficient sensitivity analysis scheme and on the 
development of good approximate models for key structural behaviors that 
maintain their accuracy for relatively large variations of values assigned to the 
design variables . 

6. If the basic studies indicate promising trends for IFM/IFMD, it will become 
important to study the program architecture of the structural analysis code that 
implements IFM/IFMD from the viewpoint of handling very large scale problems. 
I consider it important to be able to solve problems with at least 1 M degrees of 
freedom systems on a currently available computational platforms within 
reasonable resource requirements. 


I would like to congratulate Dr. Patnaik and his coworkers for bringing up this 
innovative material from abstract idea to this stage to open up a new frontier in solid 
mechanics and possible development of new class of structural analysis capabilities . 

Sincerely, 


Signed 
Hiro Miura 
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From: ”Les Berke" <lesberke§stratos . net> 

To : <Gary . R. Half ord§grc . nasa , gov> 

Subject: Patnaik* s Program 

Date: Wed, 23 Jun 1999 22:26:27 -0400 

X-Mailer: Microsoft Outlook Express 5.00.2014.211 

Dear Dr. Halford: 

I got a letter from Dr. Patnaik briefly describing his proposed long x-ange program. The 
ideas go back to the very beginning of the developments of what became the Finite 

Element Technology. I had the good fortune to be present in the trenches from the very 
begining. ' y 

Before computers structural engineers were thinking only in terms of forces, force 
equilibriums, these were, and are the quantities against which to design a structure. 
These numbers were, and are directly obtainable from analysis methods refered to as the 
Force Method. Unfortunately in case of redundant structures the equilibrium equations 
had to be augmented by a set of equations to assure compatibility of deformations under 
the obtained forces. These equations were created by a very cumbersome method of 
somewhat arbitrary and non-unique "cuts" of the structure which then were "closed" by 
the augmenting equations. J 

Because of this complexity an other method won out for computer automation, but not 
without a struggle", in the mid sixties. You can always write the same number of 
equations in terms displacements as the number of displacement unknowns in an 
articulated structure or Finite Element Model. Unfortunately forces have to be obtained 
from strains that are derivative-like quantities being obtained as small differences of 
large numbers, undercutting accuracy. The current remedy is to run inaccurate models 
to many digits of occuracy in an attempt to restore confidence in the results . 

The significant contribution of Dr. Patnaik is that in his Integrated Force Method he 
was able to develop a consistent unique set of equations based on his "global" 
compatibility conditions instead of the arbitrary artifice of "cuts" . This is a major 
and fundamental achievement that if available in the sixties, the Force Method would 
have been equally easy to automate, and would have won in that competition, forces being 
more direct and fundamental in an engineers mind. It is time to correct the past now 
that it is corractable. 

Dr, Patnaik also derived fundamental compatibility boundary conditions in continuum 
mechanics, overlooked for a Century. These conditions also play a role in the 
development of finite elements of various classes for the Integrated Force Method. 

In summary, being quite familiar with the related theories, the history of FEM 
developments and the importance and timeliness of Dr. Patnaik' s proposal I would like to 
highly recomend going ahead with his proposed developments. GLC could become a seminal 
source for a new technology in structural analysis that would eventually replace totally 
the current technology. Large investments are present in the current technology 
creating natural resistence to replace it, but there is a well known statement having to 
do with a mouse trap. , . 

Les Berke 
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NORTHWESTERN 


university 


ROBERT R. McCORMICK SCHOOL OF 
engineering and applied science 

Department of Mechanical Engineering 


August 9, 1999 

Gary R. Halford, Ph.D. 

Senior Technologist 

M.S. 49-7, NASA Glenn Research Center 
21000 Brookpark Rd. 

Cleveland, OH 44135 


Dear Dr. Halford: 


Regarding your letter of May 5, 1 found the methodology interesting; 
However, in view of the company competing needs, I cannot make any definitive 
statements about the need for IFM. 


Sincere!} 



. /ll'ist 

Ted Belytschko 

Walter P. Mun^hy Professor 
of Computational Mechanics 
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